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ABSTRACT 

Recent studies have shown that massive elliptical galaxies have total mass density 
profiles within an effective radius that can be approximated as p to t oc r~ 7 , with 
mean slope (7'} = 2.08 ± 0.03 and scatter ay = 0.16 ± 0.02. The small scatter of the 
slope (known as the bulge-halo conspiracy) is not generic in ACDM based models and 
therefore contains information about the galaxy formation process. We compute the 
distribution of 7' for ACDM-based models that reproduce the observed correlations 
between stellar mass, velocity dispersion, and effective radius of early-type galaxies in 
the SDSS. The models have a range of stellar initial mass functions (IMFs) and dark 
halo responses to galaxy formation. The observed distribution of 7' is well reproduced 
by a model with cosmologically motivated but uncontracted dark matter haloes, and 
a Salpeter-type IMF. Other models are on average ruled out by the data, even though 
they may happen in individual cases. Models with adiabatic halo contraction (and 
lighter IMFs) predict too small values of 7'. Models with halo expansion, or mass- 
follows- light predict too high values of 7'. Our study shows that the non- homologous 
structure of massive early-type galaxies can be precisely reproduced by ACDM models 
if the IMF is not universal and if mechanisms such as feedback from active galactic 
nuclei, or dynamical friction, effectively on average counterbalance the contraction of 
the halo expected as a result of baryonic cooling. 

Key words: dark matter - galaxies: elliptical and lenticular, cD - galaxies: formation 
- galaxies: fundamental parameters - gravitational lensing - stars: luminosity function, 
mass function 



1 INTRODUCTION 

Recent observations combining stellar kinematics with 
strong gravitational lensing from the Sloan ACS lens Sur- 
vey (SLACS; Bolton et al. 2006, 2008; Auger et al. 2009) 
have shown the total mass density profiles of massive ellipti- 
cal galaxies (i.e., with stellar mass above M sta r ~ 10 11 Mq or 
stellar velocity dispersion a > 200 kms" 1 ) are very close to 
isothermal (Treu & Koopmans 2004; Koopmans et al. 2006, 
2009; Auger et al. 2010b; see also, e.g., Bertin et al. 1994, 
Franx et al. 1994, Dobke & King 2006, Gavazzi et al. 2007, 
Humphrey & Buote 2010, for additional probes pointing in 
the same direction). The high precision lensing and dynami- 
cal measurements show that within the effective or half-light 
radius, the average logarithmic slope of the mass density 
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profile (ptot oc r 1 ) is (7'} = 2.078 ± 0.027 with an intrinsic 
scatter of ay = 0.16 ± 0.02 (Auger et al. 2010b). 

The physical origin of this observational result is at 
present not fully understood. We first note that the mass 
density profiles of stars and ACDM haloes are different 
from isothermal near the effective radius: If mass follows 
light (MFL), then for a de Vaucouleurs profile 7' ~ 2.3. 
In contrast, if dark matter dominates, then 7' < 1.5. The 
observational fact that 7' ~ 2 implies that both baryons 
and dark matter contribute non-negligible fractions to the 
mass within the effective radii of massive ellipticals. This is 
sometimes referred to as the "bulge-halo conspiracy" - the 
baryons and dark matter must "know" about each other so 
that the total mass profile that results from summing the 
two non-isothermal components is close to isothermal. In 
other words the observation can be qualitatively reproduced 
by combining standard ACDM halos with stellar mass den- 
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sity profiles, but only with appropriate fine tuning of the 
stellar to dark matter ratio, or equivalently of the star forma- 
tion efficiency (e.g. Keeton 2001; Gavazzi et al. 2007; Jiang 
& Kochanek 2007; Auger et al. 2010a; Nipoti, Treu & Bolton 
2008). 

Remarkably, the bulge-halo conspiracy does not happen 
at all scales. For example, the total mass density profile of 
galaxy clusters with a massive galaxy at the center is known 
to be inconsistent with isothermal (e.g., Allen et al. 2008; 
Newman et al. 2013a, b). Since dark matter halos are almost 
scale invariant in ACDM (except for mild trends in con- 
centration, e.g. Maccio et al. 2008), this difference suggests 
that understanding the origin of the bulge-halo conspiracy 
requires understanding the relevant baryonic physics. Atten- 
tion has therefore turned to cosmological hydrodynamical 
simulations of galaxy formation. 

Duffy et al. (2010) found that simulations with no 
or weak feedback can reproduce the observe 7', but they 
over predict the galaxy formation efficiencies (see also Naab 
et al. 2007). Strong feedback is needed to match the oberved 
baryon fractions, but the resulting galaxies have 7' that is 
too small (i.e., the galaxies are dark matter dominated). 
The simulations used by Duffy et al. (2010) were stopped 
at z = 2, so their conclusions are only valid if 7' does 
not evolve. Johansson et al. (2012) ran nine high resolu- 
tion "zoom- in" simulations of ~ 1O 12 M0 dark matter haloes 
down to z = (see also Remus et al. 2012). These simula- 
tions have 7' in good agreement with observations at z ~ 0, 
and also find 7' evolves from 7' ~ 3 at high redshifts to 
7' ~ 2 below z ~ 1. However, these simulations (which do 
not have strong feedback) overpredict the galaxy formation 
efficiencies by a factor of ~ 2 — 3, even assuming a Salpeter 
(1955) stellar Initial Mass Function (IMF). More recently 
Dubois et al. (2013) simulated six haloes in the mass range 
0.4— 8x 10 13 M Q both with and without feedback from Active 
Galactic Nuclei (AGN). At z = the simulations without 
AGN feedback resulted in too many baryons and overpre- 
dict the mass density slopes (7' ~ 2.3), while the simulations 
with AGN feedback reproduce the galaxy formation efficien- 
cies (assuming heavy IMFs in more massive galaxies), but 
under predict the mass density slopes (7' ~ 1.9). Thus at 
present, none of the cosmological simulations in the litera- 
ture can simultaneously match the observed 7' and galaxy 
formation efficiencies at z ~ 0. 

In this paper we aim to shed light on this issue by 
performing a quantitative comparison of flexible ACDM - 
inspired models to the observed distribution of 7' in lens 
galaxies. First, we investigate how dark matter fraction is 
required in order for models to reproduce the observed distri- 
bution of 7'. If a model has the right (i.e., observed) amount 
(and distribution of baryons) are the theoretical properties 
of ACDM haloes sufficient to explain the observed proper- 
ties of 7'? A second question we address is to what extent 
7' can help to constrain the dark matter fractions inside an 
effective radius, and hence to constrain the stellar IMF and 
dark halo response to galaxy formation. 

Observationally, a common way of measuring the dark 
matter fraction within an effective radius is achieved by 
"subtracting" the stellar mass from the total mass. Total 
masses are relatively straight forward to measure (either 
through dynamics and/or strong lensing), but stellar masses 
from spectral energy distributions are uncertain by a factor 



of ~ 3 due to the unknown form of the stellar initial mass 
function (IMF; e.g. Bell et al. 2003). 

Star counts in stellar clusters are consistent with a 
universal IMF within the Milky Way (e.g., Bastian, Cover, 
Meyer 2010). Dynamical mass-to-light ratios of extragalac- 
tic spiral galaxies are also consistent with a Milky Way IMF 
(e.g., Bell & de Jong 2001), and inconsistent with the heav- 
iei[j Salpeter (1955) IMF. Dynamical mass-to-light ratios of 
galaxies in the SAURON survey argue against a universal 
Salpeter IMF in early-type galaxies (Cappellari et al. 2006; 
Tortora et al. 2012), but they do not rule out a Salpeter-type 
IMF in all early-type galaxies. In fact, Salpeter-type IMFs 
are preferred in massive elliptical galaxies, even accounting 
for standard dark matter haloes (Treu et al. 2010; Dutton 
et al. 2011). Furthermore, there are several lines of evidence 
that require massive elliptical galaxies and even massive spi- 
ral bulges to have IMFs significantly heavier than that found 
in the Milky Way (e.g., Auger et al. 2010a; van Dokkum & 
Conroy 2010; Dutton et al. 2012; Conroy & van Dokkum 
2012; Cappellari et al. 2012; Dutton et al. 2013). 

Theoretically, the dark matter fraction within the effec- 
tive radius depends on the efficiency of galaxy formation, 
the structure of "pristine" dark matter haloes, and the re- 
sponse of the dark matter halo to galaxy formation. Galaxy 
formation efficiencies are hard to predict theoretically, but 
can be constrained through observations of weak lensing 
and satellite kinematics (e.g., Mandelbaum et al. 2006; More 
et al. 2011). In the framework of ACDM the pristine struc- 
ture of dark matter haloes is known (for a given set of cos- 
mological parameters), but the response of the dark matter 
halo is not well constrained by theory. Gas cooling or gas 
rich (dissipational) mergers are expected to make the dark 
halo contract (Blumenthal et al. 1986; Gnedin et al. 2011). 
Mass outflows (feedback) or gas poor (dissipationless) merg- 
ers are expected to make the dark halo expand (e.g., Read 
& Gilmore 2005; Governato et al. 2010; Pontzen & Gover- 
nato 2012; Maccio et al. 2012; El-Zant et al. 2001; Nipoti 
et al. 2004; Jardel & Sellwood 2009; Johansson et al. 2009). 
The relative importance of each of these processes will de- 
termine the response of the dark matter halo to galaxy for- 
mation. 

In practice, we address the two questions by using 
ACDM based models (from Dutton et al. 2012) that are 
constructed to reproduce a number of scaling relations of 
early-type galaxies including: velocity dispersion vs stellar 
mass; half-light size vs stellar mass, and dark halo mass vs 
stellar mass. The main unknowns in these models are the 
normalization of the stellar mass (and therefore the stellar 
IMF) and the response of the dark matter halo to galaxy 
formation. 

This paper is organized as follows: Our definition of 7' is 
given in §2. The constrained galaxy mass models are briefly 
discussed in §3. Results are presented in §4. The method 
used to measure ensemble average mass density slopes is 
discussed in §5. A brief discussion is given in §6, followed 



1 By heavier we mean that the IMF results in higher stellar mass- 
to-light ratios in old stellar populations. This could result from 
either an excess of low mass stars (from a bottom heavy IMF), 
or from an excess of stellar remnants (from a top heavy IMF) 
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by conclusions in §7. We adopt a ACDM cosmology with 
n A = 0.7, Q, M = 0.3 and H = 70 kms" 1 Mpc" 1 



2 DEFINITION OF MASS DENSITY SLOPE 



For an arbitrary density profile, p(r) 



- 7(r) , the local 



logarithmic slope of the density profile is dlogp/dlogr = 
— 7(r). We follow the convention of Koopmans et al. (2009) 
that 7 > 0. Thus the mass weighted slope of the density 
profile within radius r is given by 



l'(r) = 



M(r) 



-y(x)4:iTX 2 p(x) dx = 3 



Anr 3 p(r) 
M(r) ' 



(1) 



This can be expressed in terms of the local logarithmic slopes 
of the mass, M(r), and circular velocity, V(r) profiles: 



■y'(r) = 3 - 



d log M 
dlogr 



= 2 



dlogV 
dlog r 



(2) 



In the case of a power-law density profile 7'(r) = 7. But 
in general, i.e., for a non power-law density profile, 7 (r) 7^ 
7(r). Since galaxies have mass density profiles with slopes 
that decrease (i.e., become more negative) with increasing 
radius, the mass weighted slope is less than the local slope: 
7 '( r ) < l( r )- For example, for an NFW profile (which has 
an inner density slope of —1 and outer slope of —3) at the 
scale radius 7(r s ) = 2 whereas y'(r s ) ~ 1.7, and for a Hern- 
quist profile (inner slope —1, outer slope —4) at the effective 
radius "/(Re) — 2.9 whereas ~/'(R a ) ~ 2.3. In the following 
we will use the shorthand notation 7' and to indicate the 
mass weighted average slope within the effective radius. 
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Figure 1. Mean and scatter of the mass weighted density slope 
within an effective radius, 7', for models (from Dutton et al. 2012) 
and data (from SLACS — Auger et al. 2010b). The mean and 
scatter of 7' are only jointly reproduced in our model with un- 
contracted cosmologically motivated NFW haloes and a heavy 
Salpeter-type IMF (green point). Models with expansion (red 
point) or mass follows light (MFL - magenta point) over-predict 
the mean 7', while models with halo contraction (blue and black) 
under predict both the mean and scatter. 



3 MASS MODELS OF EARLY-TYPE 
GALAXIES WITH ACDM HALOES. 

This section gives a brief overview of the mass models we 
construct to reproduce the observed structural and dynam- 
ical scaling relations of early-type galaxies in the Sloan 
Digital Sky Survey (SDSS; York et al. 2000; Abazajian 
et al. 2009). A more detailed discussion is given by Dutton 
et al. (2012). Our mass models consist of spherical distribu- 
tions of stars and dark matter. The stars are modelled as the 
sum of two Sersic profiles (with n = 1 and n = 4), while the 
dark matter is an NFW (Navarro, Frenk, & White 1997) pro- 
file with "pristine" concentration parameters from Maccio 
et al. (2008) and modified by halo contraction/expansion 
following Dutton et al. (2007). 

The stellar profiles of the models are constructed 
to reproduce the size-mass relation of early-type galaxies 
from Dutton et al. (2012) which uses sizes from Simard 
et al. (2011), and stellar masses from the MPA/JHU grouj0 
which assume a Chabrier (2003) IMF. The halo masses 
are chosen to follow the observed relation between stel- 
lar mass and halo mass for early-type galaxies from Dut- 
ton et al. (2010) which combines results from weak lensing 
(e.g., Schulz et al. 2010) and satellite kinematics (e.g., More 
et al. 2011). The models have three sources of scatter: scatter 
in galaxy size at fixed stellar mass (constrained by observa- 
tions: Dutton et al. 2012); scatter in dark halo mass at fixed 

2 Available at http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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stellar mass (constrained by observations: More et al. 2011); 
and scatter in dark halo concentration at fixed dark halo 
mass (constrained by theory: Maccio et al. 2008). 

The velocity dispersion vs stellar mass (a.k.a. Faber 
& Jackson 1976) relation is used to constrain the allowed 
combinations of stellar IMF and dark halo response in 
these models. We consider five different halo responses 
ranging from standard adiabatic contraction (Blumenthal 
et al. 1986) to unmodified NFW haloes to mass-follows-light 
(i.e., maximum expansion). In order to reproduce the slope 
and zero point of Faber- Jackson relation all these models 
require "heavier" IMFs in more massive galaxies, and (triv- 
ially) "heavier" IMFs in models with stronger halo expan- 
sion. These set of models are physically realistic (at least in 
terms of reproducing the scaling relations) while having a 
wide range of dark matter fractions and inner dark matter 
halo density profiles. 

The correlation between the scatter in the velocity-mass 
relation with the scatter on the size-mass relation — equiv- 
alent to the tilt of the Fundamental Plane (FP; Dressier 
et al. 1987; Djorgovski & Davis 1987; Ciotti et al. 1996) - 
was used by Dutton et al. (2012) as an additional constraint 
to distinguish between these models. For galaxies in the mass 
range 10 10 < Msps ^ 10 n MQ, only models in which mass- 
follows-light were able to match this additional constraint. 
For the most massive galaxies (and those which are relevant 
for this paper) models with unmodified NFW haloes were 
favored over MFL and adiabatic contraction models. 
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Figure 2. Distribution of mass weighted density slope within an effective radius, 7', for our models compared with the observed (intrinsic 
distribution) from Auger et al. (2010b, black dotted line). The histograms show 7' for the total mass (black, solid), stellar mass (blue, 
short-dashed), and dark matter (red, long-dashed). Each panel has a model with a different dark halo response: Adiabatic contraction, 
v = 1.0 (upper left); mild contraction v = 0.5 (upper right); no contraction v = 0.0 (lower left); and expansion v = — 0.5( lower right). 
The model with mass follows light is equivalent to that of the stellar component. Only the model with uncontracted cosmologically 
motivated dark matter haloes and a heavy IMF reproduces the observations. 



4 RESULTS 

In order to make a fairer comparison between our mod- 
els and data from SLACS we consider model galaxies that 
match the distribution of velocity dispersions for the SLACS 
lenses. Specifically, we select model galaxies with a log- 
normal filter on velocity dispersion: mean log(CT/[kms -1 ]) = 
2.40 and standard deviation 0.07. With this selection 
the average offset between the model stellar masses and 
those obtained assuming a Chabrier IMF are Aimf = 
log 10 (M st ar/M S ps) = (0.08, 0.20, 0.27, 0.31, 0.33) for halo re- 
sponses of v = (1.0,0.5,0.0, -0.5, MFL), respectively 0. For 
reference a Salpeter IMF corresponds to Aimf — 0.23, and 
thus models with uncontracted dark matter haloes (u = 0.0) 
and expansion [y — —0.5) have IMFs heavier than Salpeter. 

The observed 7' are not strictly measured within any 
uniform aperture such as one effective radius, Re- The aper- 
ture depends on the Einstein radius and the physical aper- 

3 For easy comparison with previous studies Aimf = log a, 
where a is the so-called IMF mismatch parameter, e.g., Treu 
et al. (2010), Dutton et al. (2013). 



ture of the SDSS fiber. The former depends on the velocity 
dispersion of the lens and the redshifts of the lens and source, 
while the latter just depends on the redshift of the lens. How- 
ever, for SLACS lenses the aperture is typically between one 
and one-half an effective radii (e.g., Koopmans et al. 2009). 
To test whether 7' depends on the radius within which it 
is measured we have compared 7' within R c and R e /2. In 
general y'(R c /2) < "/'(R a ), but the changes are small. The 
changes in the median 7' for our five models are as fol- 
lows: j'(Re) -7' (.Re/2) = (0.028,0.022,0.001,0.070,0.141) 
for halo responses of v = (1.0, 0.5, 0.0, -0.5, MFL). Thus the 
mass density slopes varies between one and one-half an effec- 
tive radius by much less than the intrinsic scatter, except for 
the MFL model, which we will show below is a poor match 
to the data for other reasons (and is well known to be a poor 
match to massive elliptical galaxies anyway) . For simplicity, 
in what follows we measure the mass density slopes within 
one effective radius, with no significant loss of precision. 

Our main results are summarized in Fig. \T\ which shows 
the scatter in 7' vs the mean of 7' for our 5 models and 
observations from Auger et al. (2010b). The distributions of 
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7' for the four models with dark matter haloes are shown in 
Fig. [2] The model with MFL is only comprised of stars and 
therefore has 7' = 2.33^ 



1-0.02 
-0.03- 



2.5 



4.1 Average mass density slopes 

All of the models broadly reproduce the close to isothermal 
density profiles that are observed, with average slopes in the 
range f .8 < 7' < 2.4. However the model with uncontracted 
NFW haloes (and a slightly heavier than Salpeter IMF) pro- 
vides the best match to the observed average 7' (Fig. [T] and 
yellow shaded region in Fig. [5J). This is in good agreement 
with the conclusions of Auger et al. (2010a) who favored 
uncontracted haloes over any type of halo contraction (they 
did not consider expansion or MFL), and IMFs significantly 
heavier than Chabrier. It is also in excellent agreement with 
the, completely independent, constraints from the tilt of the 
fundamental plane from Dutton et al. (2012). 



4.2 Scatter in mass density slopes 

Most importantly however, our models also reproduce the 
observed scatter of the distribution, i.e. the tightness of the 
bulge-halo conspiracy. The distribution of 7' in the models 
is roughly Gaussian (although slightly skewed) with scatter 
smaller or equal to the observed intrinsic scatter of 0.16 ± 
0.02. As with the mean 7', the model with uncontracted 
NFW haloes best matches the scatter with 7' = 2.06^° 

Given that there is significant overlap between the ob- 
served distribution of 7' and all of our models, it is entirely 
possible that the full range of halo responses occurs in real 
galaxies. Indeed, there are already observational hints that 
this is the case. While on average a Salpeter-type IMF fa- 
vors uncontracted NFW haloes (e.g., Treu et al. 2010; Auger 
et al. 2010a; Dutton et al. 2011, Dutton et al. 2012), the de- 
tailed study of a massive elliptical, enabled by the presence 
of a double Einstein Ring, favors a Salpeter-type IMF and 
a contracted halo (Sonnenfeld et al. 2012). From a theoreti- 
cal standpoint this could be connected to the specific merger 
history of each galaxy (e.g., Nipoti et al. 2012 and references 
therein) or perhaps to the mode of starformation and re- 
sulting IMF (Hopkins 2012). Below we will explore whether 
7' correlates with other observables. This could potentially 
provide some clues as to the origin of each galaxy. 



4.3 Mass density slopes of stars and dark matter 

In Fig. [2] the red histograms show the mass weighted den- 
sity slopes within R c of the dark matter haloes, 7dm- These 
range from ~ 1.7 for adiabatic contraction, to ~ 0.8 for 
expanded haloes. The blue histogram shows the the mass 
weighted density slopes within R c of the stars, these have 
7 ^ tar ~ 2.33 ± 0.02. 

The relative contribution of the stars and dark matter 
to the total 7' is given by the dark matter fraction within 
an effective radius, /dm (-Re). As expected, for a given halo 
response model 7' is strongly anti-correlated with the dark 
matter fraction (Fig. [3J. The lines show the median 7' as a 
function of dark matter fraction, while the dots show the me- 
dian 7' and /dm (Re). While for an arbitrary halo response 
7' does not uniquely predict the dark matter fraction, it 
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Figure 3. Correlation between mass density slope, 7', and dark 
matter fraction within an effective radius, /dm, for our models. 
The colored lines show the median mass density slope as a func- 
tion of the (spherical) dark matter fraction, while the dots show 
the median values. For a given halo response the mass density 
slope is determined by the dark matter fraction. 



does place some useful limits. For example, if 7' > 2, then 

/DM (-Re) < 0.2. 

The average dark matter fractions vary from ~ 0.55 for 
adiabatic contracted NFW haloes (black dot) to ~ 0.10 for 
expanded NFW haloes (red dot). For NFW haloes (green 
dot), the average dark matter fraction is just ~ 20%, but 
this is enough to lower the mean 7' by 0.27 and increase the 
scatter in 7' by 0.14 relative to the MFL model. 



4.4 What does 7' correlate with? 

Previous observational studies have shown that 7' is largely 
uncorrelated with many physical properties of galaxies, and 
does not evolve significantly since z ~ 1 (e.g., Koopmans 
et al. 2009; Auger et al. 2010b; see however Ruff et al. 2011 
and Bolton et al. 2012 for tentative detections of a mild 
evolutionary trend). The most significant correlation found 
by Auger et al. (2010b) was between 7' and total surface 
density, with denser galaxies having higher 7' . Such a corre- 
lation is expected because denser galaxies have either higher 
density stellar components (which increases the contribution 
of the stars to 7') and/or stronger halo contraction (which 
results in higher 7dm)- This corelation might provide an im- 
portant clue to understanding the origin of the small vari- 
ations from galaxy to galaxy. The central density could be 
a by product of the epoch of formation of the inital core 
and successive merger history and perhaps be related to the 
conditions of the gas that created the bulk of the old stars. 

Fig. [4] shows correlations between 7' and a number of 
physical galaxy parameters: stellar surface mass density; ef- 
fective radius; stellar mass; and stellar velocity dispersion. 
Lines show the median relations for the models, the points 
show the model galaxies with uncontracted NFW haloes, 
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Figure 4. Correlations between mass density slope within the effective radius, 7', and galaxy observables: stellar surface density (£ e , 
upper left), effective radius (i? c , upper right), stellar mass (Mgps, lower left), and stellar velocity dispersion (<j e 2, lower right). The models 
are given by colored lines: adiabatic contraction (v = 1.0, black long-dashed); mild contraction (y = 0.5, blue dotted); no contraction 
{v = 0.0, green solid); halo expansion (1/ = —0.5, red short-dashed). In addition, the points show individual galaxies for the no contraction 
model (which provides the best match the the observed mean and scatter of 7' - see Fig.[TJ. The scatter about the median relation for 
this model, cry , is given in the lower right corner of each panel. The tightest correlation in the models is between 7' and stellar density, 
while the weakest is with velocity dispersion. The data from SLACS (Auger et al. 2010b) are given by filled symbols with error bars. For 
both models and data the stellar masses and surface densities have been scaled to a Salpeter IMF. The trends seen in the no contraction 
model are also seen in the observations, providing further support for its validity. 



and the filled points with error bars show the observations 
from SLACS (Auger et al. 2010b). To ease comparisons be- 
tween the different models and the data the stellar masses 
and stellar densities have been re-scaled to a Salpeter IMF. 
The most significant correlation in the data is with stellar 
density, followed by effective radius and stellar mass. These 
trends are qualitatively reproduced by the models. To the 
eye, the data appear to suggest a weak positive correlation 
with velocity dispersion, which is of opposite sign to the 
models. However, the measurement errors on 7' are posi- 
tively correlated with the errors on rr C 2 (Auger et al. 2010b), 
and therefore the trend is statistically insignificant (slope 
0.07 ± 0.08) and consistent with the weak opposite trend 
seen in our models. 

Thus, not only does our model with uncontracted NFW 
haloes reproduce the mean and scatter of the observed 7', it 
also reproduces the correlations between 7' and other galaxy 
observables. We note this achieved without any arbitrary 



fine-tuning of the distribution of baryons and dark matter 
in galaxies. 

In our models, for a given halo response, there is almost 
no scatter between 7' and dark matter fraction. Thus the 
strength of the correlation between 7' and various galaxy 
parameters are largely determined by the strength of the 
correlations with dark matter fraction (see Fig. [5j. 



5 ENSEMBLE MASS DENSITY SLOPE 

In this section we comment on a method to infer the average 
mass density slope of the dark matter within the effective 
radii of massive early-type galaxies. This method was re- 
cently used by Grillo (2012) to argue in favor of very steep 
dark matter density profiles, which is in apparent contradic- 
tion to our results (as well as those by Auger et al. 2010a). 

The main observables are the total projected mass, 
A/tot (Jfein), within the Einstein radius, ifein, obtained from 
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Figure 5. Correlations between dark matter fraction within the effective radius, /dm(Kc), and galaxy observables (as in Fig. [4} for our 
models. The scatter about the median relation for the no contraction model, try , is given in the lower right corner of each panel. The 
size of the scatter in these relations is directly correlated with the size of the scatter in the equivalent relations in Fig. [4] 



strong gravitational lensing. For a given IMF one can cal- 
culate the total mass in stars, Msps, and the mass in stars 
within the Einstein radius, Msps(< ifein). Thus for each 
lens, assuming an IMF, one can calculate the dark matter 
surface density within the Einstein radius: 

SDM^Em) = [Mtot(< #Ein) - M SPS (< #Ein)] /^lin- (3) 

Since the Einstein radii are in general different for each lens 
system, one can determine the average dark matter den- 
sity profile for an ensemble of gravitational lenses (see, e.g., 
Rusin & Kochanek 2005 for a similar strategy). However, 
the lenses have different masses and sizes, and therefore the 
esemble average cannot be conducted in physical scales. A 
useful approach is to scale the size and density measure- 
ments with the intrinsic scales of the galaxies: the effective 
radii R e , and effective surface densities E c = Msps/R^. In 
this way, one can then define dimensionless sizes and dark 
matter densities via 

A = Rvin/Re, (4) 

and 

* = EDM(i?Ei n )/S c (i? c ). (5) 

Plotting $ vs A thus yields an ensemble density profile of 



the dark matter. The correlation can be described with a 
power-law: 

* = ahP (6) 

to yield an effective de-projected dark matter density slope 
7DM,g = -/3 + 1. 

Applying this method to 39 massive early-type lenses 
from the SLACS survey, Grillo (2012) found /3 = -1.04±°;^ 
(corresponding to 7dm, g = 2.04lg' jg) for a Chabrier IMF, 
and p = -0. 771^7 (corresponding to 7dm. g = l-77to.62) 
for a Salpeter IMF. Since an uncontracted NFW halo will 
have P ~ —0.1 (7dm — 1-1), this provides marginal evidence 
for halo contraction in response to galaxy formation (assum- 
ing the IMF is lighter than Salpeter), in contradiction to our 
results. However, as we discuss below, there are two issues 
that explain this apparently discrepancy. 

The first issue is that the IMF is likely Salpeter or 
heavier, which will obviously require shallower dark mat- 
ter slopes. For example Auger et al. (2010a), Dutton 
et al. (2012), and Conroy & van Dokkum (2012) all favor 
IMFs Salpeter or heavier in the most massive galaxies. We 
note that for a Salpeter IMF or heavier our results are in 
agreement with Grillo's. 

The second issue is more subtle and has to do with an 
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implicit assumption of the method used by Grillo (2012). 
Grillo assumes that the re-scaling of the dark matter den- 
sities and Einstein radii does not introduce spurious cor- 
relations. Specifically, in order for ^(A) to be equivalent 
to £DM(i?Ein) requires that Esps/Rf is independent of 
R E in- Using SLACS data from Auger et al. (2010b) we find 
that this is not the case. Fitting a relation of the form 
E e /-Re K (-REin) 5 on a grid of f3 over the interval [—2.5, 1], 
and then fitting a linear relation between <5 and /3 we find 
that 5 = -0.078-0.586,3. Thus the true dark matter density 
slope is given by 

7dm = 7dm, g - * = 0.664 + 0.414 7DM , g . (7) 

So a measured isothermal mass slope (i.e., 7dm, g — 2) im- 
plies a shallower true mass slope of 7dm — 1-49. 

Correcting the results from Grillo (2012) using Eq. [7] 
yields: 7 fjM.g — l-5liolii for a Chabrier IMF, and 7 fjM, g = 
1.40^0 26 fo r a Salpeter IMF. In the context of our mod- 
els, this would favor weak contraction (y ~ 0.5), which is 
similar to that found in cosmological simulations by Abadi 
et al. (2010). Allowing for an IMF slightly heavier than 
Salpeter (as we favor) would increase the inferred /3 to be 
in even closer agreement with our results. In summary, we 
find that there is no conflict between the ensemble average 
dark matter density slope as derived by Grillo (2012) and 
our results based on the fundamental plane and total mass 
density slopes 



6 DISCUSSION 

We have shown that the observed distribution of mass den- 
sity slopes can be reproduced, precisely, if galaxies with re- 
alistic sizes and concentrations are embedded in cosmolog- 
ically motivated dark matter haloes. We note that a key 
requirement for our models to reproduce the observed 7' is 
that the concentration of the stars needs to be high enough 
so that 7star(-Rc) > 2 — because ACDM haloes on their own 
never have such steep density profiles on such scales. For 
example, galaxies with exponential surface density profiles 
(i.e., characteristic of galaxy disks) would be unable to result 
in 7' (R c ) ~ 2. 

Our model assumes that the scatter in galaxy sizes at 
fixed stellar mass is uncorrelated with the structure of the 
dark matter halo. This assumption was made for simplic- 
ity (and remarkably it seems to work), but it should be 
tested using galaxy formation models. For example, for spi- 
ral galaxies we do expect correlations between the scatter in 
galaxy sizes with scatter in halo concentrations and galaxy 
formation efficiencies (see Fig. 7 in Dutton et al. 2007). It re- 
mains to be seen whether such correlations exist for spheroid 
dominated galaxies. 

Our results favor unmodified cosmologically motivated 
(Navarro et al. 1997; Maccio et al. 2008) dark matter haloes. 
This is somewhat surprising, since there are many processes 
that should modify the dark halo structure (e.g., gas ac- 
cretion, gas outflows, minor and major mergers). One pos- 
sibility is these various processes are occurring stochasti- 
cally, and cancel each other out on average. Alternatively, 
we might expect that gas accretion dominates at early times, 
followed by gas outflows, then dry mergers. Under this sce- 
nario, we would expect the halo to initially contract, and 



then to slowly undo this contraction over time. Thus a mea- 
surement of the evolution of the halo response would be able 
to distinguish between these formation scenarios for massive 
elliptical galaxies. This should be possible with the current 
generation of galaxy scale lenses extending to higher red- 
shifts than SLACS (Sonnenfeld et al. 2013, in preparation).. 



7 CONCLUSIONS 

In the past few years observations have shown that massive 
early-type galaxies have close to isothermal density profiles 
within one effective radius (Koopmans et al. 2009; Auger 
et al. 2010b). We use ACDM based mass models constructed 
to reproduce a number of the observed scaling relations of 
early-type galaxies (including the Faber-Jackson relation) 
to address two questions: 1) For models with the correct 
distribution of baryons, do the properties of ACDM haloes 
result in models that reproduce the observed properties of 
7'? 2) Do the observed properties of 7' help to distinguish 
between models with different stellar Initial Mass Functions? 
We summarize our results as follows: 

• The median values of the mass density slopes of all our 
models are roughly isothermal: 1.8 < 7' < 2.4. The observed 
mean value of (7') = 2.08 ± 0.03 is best reproduced by a 
model with uncontracted NFW haloes and stellar masses 
~ 0.27 dex higher than that obtained assuming a Chabrier 
(2003) IMF (Fig. P. 

• The scatter in 7' in our models with dark matter haloes 
is small 0.10 < ay < 0.17. The observed (intrinsic) scatter 
of ay = 0.16 ± 0.02 is best reproduced by our model with 
uncontracted dark matter haloes (Fig. [1} . 

• In observations the mass density slope correlates with 
stellar surface density and effective radius, with higher 7' in 
higher density and smaller galaxies. As with other properties 
of 7' this correlation is best reproduced by our model with 
uncontracted NFW haloes (Fig. 2]). 

• In the models the tightness of the correlations involving 
7' are largely determined by the tightness of the correspond- 
ing correlations involving (spherical) dark matter fraction 
within the effective radius (Fig. [SJ . 

In conclusion, our study shows that the many of the 
observed properties - including the bulge-halo conspiracy 
and other classic scaling relations like the Faber Jackson 
and size mass relation - of massive early-type galaxies can 
be precisely reproduced by ACDM models under two con- 
ditions: i) the IMF is not universal; ii) mechanisms such as 
feedback from active galactic nuclei, or dynamical friction, 
effectively counterbalance on average the contraction of the 
halo expected as a result of baryonic cooling. We empha- 
size that no correlations between the scatter in galaxy and 
dark halo parameters (or equivalently, no fine-tuning be- 
tween the baryons and dark matter) is needed in order for 
our models to reproduce the observed small scatter in total 
mass density profiles. Although our models are clearly an 
oversimplified description of reality, they demonstrate that 
self-consistent ACDM -inspired models can be found. This 
gives hope that the recent improvements in numerical sim- 
ulations including baryonic physics could soon lead to to 
realistic models for the formation of massive type galaxies 
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in quantitative agreement with the tight constraints pro- 
vided by observations. In turn this might help explain what 
is the relative contribution of the different processes that 
contribute to produce the effective profiles that we observe 
Our work also provides further evidence against a universal 
IMF for galaxies. In agreement with earlier studies (Auger 
et al. 2010a; Spiniello et al. 2011; Sonnenfeld et al. 2012, 
van Dokkum & Conroy 2010; Conroy & van Dokkum 2012; 
Smith et al. 2012; Spiniello et al. 2012; Dutton et al. 2012; 
Cappellari et al. 2012; Dutton et al. . 2013) we find that a 
"heavy" Salpeter-like IMF is preferred for massive galaxies 
over the light Chabrier-like IMFs usually preferred for Milky 
Way type galaxies. The origin of this non-universality is still 
hotly debated. However, the fact that the only quantity that 
seems to correlate with the mass structure of massive early- 
type galaxies is the central surface density could be a clue 
that the density of stars (and therefore possibly gas at the 
epoch of peak star formation) is an important ingredient 
of shaping the final stellar initial mass function (Hopkins 
2012). 
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